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ABSTRACT 

We report on deep near-infrared F125W (J) and F160W (H) Hubble Space Telescope Wide Field 
Camera 3 images of the z = 6A2 quasar J1148+5251 to attempt to detect rest-frame near-ultraviolet 
emission from the host galaxy. These observations included contemporaneous observations of a nearby 
star of similar near-infrared colors to measure temporal variations in the telescope and instrument 
point spread function (PSF). We subtract the quasar point source using both this direct PSF and a 
model PSF. 

Using direct subtraction, we measure an upper limit for the quasar host galaxy of mj > 22.8, 
ttih > 23.0 AB mag (2 a). After subtracting our best model PSF, we measure a limiting surface 
brightness from 0'/3 — 0'/5 radius of \ij > 23.5, jin > 23.7 AB mag arcsec -2 (2 a). We test the ability 
of the model subtraction method to recover the host galaxy flux by simulating host galaxies with 
varying integrated magnitude, effective radius, and Sersic index, and conducting the same analysis. 
These models indicate that the surface brightness limit (fij > 23.5 AB mag arcsec -2 ) corresponds to 
an integrated upper limit of mj > 22 — 23 AB mag, consistent with the direct subtraction method. 
Combined with existing far-infrared observations, this gives an infrared excess \og(IRX) > 1.0 and 
corresponding ultraviolet spectral slope j3 > —1.2 ±0.2. These values match those of most local 
luminous infrared galaxies, but are redder than those of almost all local star- forming galaxies and 
Lyman break galaxies. 

Subject headings: galaxies: high-redshift, methods: observational, quasars: individual (SDSS 
J1148+5251) 



1. INTRODUCTION 

Since their discovery bv iSchmidtl (|1963f ). quasars have 
been the best and most easily observed beacons to probe 
the distant universe. The unification mode l for active 
galactic nuclei (AGN, e.g., lAntonuccl I1993L and refer- 
ences therein) views quasars as unobscured active nuclei, 
with a super-massive black hole (SMBH) as the central 
power-house behind the dominan t nonstellar cont i nuum . 

Ground-based studie s (e.g., iBoroson fc Okd [1982, 
iMcLeod fc Riekel Il994f ) are limited by atmospheric 
seeing, and show that observing the underlying stellar 

1 School of Earth and Space Exploration, Arizona State 
University, P.O. Box 871404, Tempe, AZ 85287, USA 

2 Department of Physics, Arizona State University, P.O. Box 
871504, Tempe, AZ 85287, USA 

3 Space Telescope Science Institute, Baltimore, MD 21218, 
USA 

4 Steward Observatory, The University of Arizona, Tucson, 
AZ 85721, USA 

5 Carnegie Observatories, 813 Santa Barbara Street, 
Pasadena, CA 91101, USA 

6 Department of Physics and Astronomy, University of 
Alabama, Box 870324, Tuscaloosa, AL 35487, USA 

7 Leiden Observatory, Leiden University, PO Box 9513, 2300 
RA Leiden, The Netherlands 

8 Department of Astronomy and Astrophysics, The Pennsyl- 
vania State University, 525 Davey Laboratory, University Park, 
PA 16802, USA 

9 Institute for Gravitation and The Cosmos, The Pennsylva- 
nia State University, University Park, PA 16802 

10 p r i nce ton University Observatory, Princeton, NJ 08544, 
USA 

11 Department of Physics and Astronomy, The University of 
Missouri, 701 South College Ave, Columbia, MO 65211, USA 



popul ations is not trivial. iTargett. Dunlop. fc McLurd 
(2012) have recently imaged z ~ 4 quasar host galaxies 
using adaptive optics with the European Southern 
Observatory Very Large Telescope, but these required 
the best seeing conditions and targets with nearby 
(< 40 arcsec) bright stars. Studies with the Hubble 
Space Telescope (HST) have concentrated on quasars at 
z ~ 0.5 — 3, yielding important constraints on the mor- 
phology, lum inosity, and stellar populations of quasar 
hosts (e.g., iDisnev et all H995j [M cLeod fc Riekd Tl995l 
Bah call et all I1997L iMcLure et al. 1999, Ridg wav et atl 
I2001L iH utching s et al.l I2002L [Dunlop et aLM 2Q03. 
Peng et al. 2006, Zakamska et al. 2006). 

iGoto et all (j2009f ) and iWillott et all ([201 lh used 
ground-based telescopes to detect extended Lyman-a 
emission out to 15 kpc around the z = 6.42 quasar 
CFHQS J232908-030158. This underscores the need 
to observe at longer, rest-frame near-ultraviolet wave- 
lengths if emission from the host galaxy is to be differ- 
entiated from surrounding gas. 

Imaging quasar host galaxies allows one to examine 
how mergers, starbursts, and other changes in environ- 
ment affect the central nucleus. Observations of molecu- 
lar gas and the far-infrared continua of z ~ 6 quasars have 
suggested extrem ely high average s tar formation rates of 
> 1000 M© yr _1 ([Wang et al.ll2010f ). implying potentially 
luminous hosts. As many quasars are hosted in massive 
galaxies, imaging z~6 quasar hosts may reveal the prop- 
erties of the most massive first galaxies. The enormous 
SMBH mass (> 10 9 M Q ) associated with Sloan Digital 
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Sky Survey (SDSS, lYork et all I2000D quasars at z ~ 6 
represents a great theoretical challenge that heavily con- 
strains possible formation methods an d merger histories 
(jVolonteri fe Reesll2QQ6l iLi et al.ll2QQ7h . 

SDSS J114816.64+525150.3 (hereafter J1148+5251), is 
the best-studied member of the z~6 quasar population, 
having been extensivel y observed a t mult iple wavelengths 
since its discov ery bv iFan et al.l (I2QQ3I) . Near- infrared 
spectr oscopy by IWillott et al.l (j2003[ ) and llwamuro et al.l 
(|2004l ) measured the Mg II and Fe II features, estimat- 
ing a mass of 3 x 10 9 M for the SMBH, an accre- 
tion rate near the Eddington limit, and an Fe Il/Mg II 
ratio consistent with quasars at lo wer redshifts. Ra- 
dio observations of CO lines (e.g., I Walter et al.l 2003, 
iRiechers et al.ll2009[ ) indicate the presence of 2.2 — 2.4 x 
10 10 Mq of high-excitation molecular gas extending to 
~ 2.5 k pc (r = f !4b) Stu d ies of the fC III l ine at 
158/xra (jMaiolino et al.l I2005L iWalter et all I2009D pro- 
vide evidence that the quasar host galaxy is undergo- 
ing a vigorous starburst, with an estimated star forma- 
tion rate density of ~ 1000 M yr _1 kpc -2 extending 
over kiloparsec scales. Studies of the continuum emis- 
sion at far-infrared (FIR) wavelengths indicate a warm 
dust component with an AGN-correct ed FIR luminos- 
ity of 9.2 x 10 12 L Q (Wa ng et~aTll2010f ) a nd correspond- 
ing dust mass of 4 2 - 7.0 x 10 8 dBertoldi et al.l 
I2QQ3L iRobson et all 12004 iBeelen et alll2QQ6h . Locally, 
this is in the range of ultra-luminous infrared galax- 
ies (ULIRGs, galaxies with Lpm > 10 12 £©). Near- 
an d mid-infrared Sp itzer Space Telescope observations 
by Uiang et ah! (j2QQ6h show clear evidence for prominent 
hot dust within the galaxy. All of these observations 
argue for a host galaxy with a significant stellar mass 
component undergoing an extreme episode of star for- 
mation. 

In this Letter, we report on near-infrared imaging of 
J1148+5251 with the HST Wide Field Camera 3 (WFC3) 
and our methods for characterizing and subtracting the 
instrument and telescope Point Spread Function (PSF). 
In $2j we describe our near-infrared WFC3 imaging of 
Jl 148+5251 and a nearby star, used for PSF character- 
ization. In $3j we describe our method for subtracting 
the central point source from the quasar images. We 
describe our simulations for assessing the reliability of 
our subtraction method in $4] Finally, in $5j we dis- 
cuss the implications of these results and our plans for 
future investigation. For this paper, we adopt a ACDM 
cosmology with Hp = 70.3 km s" 1 Mp c -1 , Qm = 0.271, 
and = 0.729 (jKomatsu et al.ll2011h . U nless other wise 
stated, all magnitudes use the AB system (Oke 1974) and 
h aye been corrected for Galactic extinction using the map 
of lSchlegel etall (fl998h . 

2. OBSERVATIONS AND DATA REDUCTION 

The HST observations of Jl 148+5251 were performed 
on 2011 January 31 (HST Program ID 12332, PI: R. 
Windhorst) using the WFC3 IR channel with the F125W 
(Wide J) a nd F160W (WFC3 H) filters. Previous pro- 
grams (e.g.. lHutchings et al.ll2002f ) found that the quality 
of empirical quasar point source subtractions was signif- 
icantly affected by PSF time variability. This variability 
is mostly due to so-called "spacecraft breathing" effects, 
i.e., thermally induced defocus of the HST secondary mir- 
ror due to movement of the Optical Telescope Assembly 



TABLE 1 

Summary of Observations 



Target 


Filter 


Exposure 
Time (s) 


S/N 


SDSS J1148+5251 


F125W 


2478 


2400 


SDSS J1148+5251 


F160W 


3646 


3760 


2MASS J11552259+4937342 


F125W 


208 


1730 


2MASS J11552259+4937342 


F160W 


335 


2200 



as the telesco pe goes into and out of Earth shadow (e.g., 
iHers hev 1998). Examining these studies and instrument 
reports, we identified two primary sources of thermal 
variation that we expected to affect our observations - 
the spacecraft attitude and the orbital day-night cycle. 

To minimize thermal variations due to spacecraft atti- 
tude, we constrained our PSF star to be within 5° of the 
target quasar, and required that the quasar and PSF star 
observations be collected in contiguous orbits. We fur- 
ther requested (and were granted) that our observations 
be scheduled immediately following another observation 
near the same celestial coordinates, thus minimizing 
thermal equilibration time in our first orbit. We selected 
the PSF star to match the quasar near-infrared colors 
((J-H) = 0.55 Vega mag, (H-K') = 0.72 Vega mag) as 
closely as possible, to minimize differences in wavelength- 
dependent PSF features. Many stars with colors similar 
to the quasar also have SDSS spectra, since they were tar- 
geted by SDSS as high-redshift quasar candidates. We 
examined these spectra where available, to reject obvi- 
ous spectroscopic binaries or other contaminants. Af- 
ter checking the resulting candidate stars for HST guide 
stars, we selected the star 2MASS Jl 1552259+4937342 
(spectral type K7, mj = 16.272 ± 0.079 Vega mag, 
(J-H) = 0.488 ± 0.158 Vega mag, (H -Ks) = 0.651 ± 
0.176 Vega mag. iSkrutskie et aLll2QQ6h as our final PSF 
target. 

Unfortunately, there is no way to eliminate thermal 
variation due to the orbital day /night cycle. Thus, we 
constructed our orbits to ensure that we fully sampled 
this cycle, and that equal fractions of the final, com- 
bined quasar and PSF star images would come from 
any given location in orbital phase. Each quasar ex- 
posure was matched by several shorter PSF star expo- 
sures, taken at the same sub-pixel dither point at the 
same phase within the orbit. This phase-space sampling 
is summarized in the bottom panel of Figure [U Any 
given detector location was never exposed beyond half- 
well depth in a single orbit, thus avoiding saturation or 
detector persistence. The quasar was observed in this 
pattern for three orbits in total, and the PSF star was 
observed for a single orbitP^l. The total exposure times 
for each object and filter combination are summarized 
in Table [U Anal ysis was performed on Multi drizzle- 
combined images ( Koekemoer et aLll2002l l2011f ) with an 
output pixel scale of 0'/065, to achieve Nyquist sampling 
of the PSF core in both filters, enabling accurate spatial 
shifting. The cosmic ray rejection step of Multidrizzle 
was disabled, since the WFC3 IR MULTIACCUM readout 

12 An additional PSF star, 2MASS J11403198+5620582, was 
also observed for a single orbit to allow for inter-orbit interpolation 
of the PSF measurement. This observation suffered from a poor 
guide star acquisition and was unusable. 
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Fig. 1. Empirical PSF subtraction. Left panels: Multidrizzle-combined F125W (J, top) and F160W (if, bottom) WFC3 images 
of Jl 148+5251 after sky subtraction. Pixels are 0^065. Middle Panels: Scaled and shifted PSF star, as fit by GalFit. Right panels: 
Fit residuals, showing a net positive residual flux, but high noise. We measured the integrated residual flux using a 0"5 radius aperture 
(gray circle), obtaining upper limits of mj > 22.8 mag, rnn > 23.0 mag (2 a). All images are displayed with the same logarithmic stretch. 
Bottom panels: PSF star and one quasar orbit, highlighting the relative phasing of corresponding dither points to compensate for spacecraft 
breathing. Exact phase matches are not possible due to buffer dumps and specific readout sequences. Corresponding dither points were 
centered at similar positions on the detector to account for field-dependent PSF variability. 

mode provided sufficient cosmic ray rejection. 



3. POINT SOURCE SUBTRACTION 

To subtract t he quasar point s ource, we used the pro- 
gram GalFit (jPeng et all I2Q1QD to fit a PSF single- 
component model to the quasar image. We transformed 
the Multidrizzl e-generated weight ma ps into uncer- 
tainty maps as in iDickinson et all ([20041 ), including the 
effects of correlated noise and shot noise from the quasar, 
and supplied these to GalFit as the pixel-to-pixel uncer- 
tainty ("sigma") image. We then subtracted the best-fit 
model from the original image, and inspected the resid- 
ual. 

We first attempted this subtraction using the image of 
our PSF star as the model PSF. The results are shown 
in Figure [TJ We measured the residuals using a 0"5 ra- 
dius aperture, obtaining upper limits of mj > 22.8 mag, 
rriH > 23.0 mag (2 a). This includes the total noise con- 
tribution from both the quasar and the empirical PSF, 
measured by scaling the PSF uncertainty map by the 
same factor as in the fit, and adding it in quadrature 
to the quasar uncertainty map. The noise contribution 
from the subtracted PSF is comparable to that of the 
quasar since the two images have comparable S/N (see 
Table [TJ, which leads to cosmetic defects (holes) in the 
subtraction, despite the net positive residual. 



We also generated a TinyTin0 (|Krist, Hook, fc Stoehrl 
2011) model of the PSF, which we calibrated to the 
the PSF star observations. A 5x spatially oversam- 
pled TinyTim model was generated for each WFC3 ex- 
posure to allow for subpixel shifti ng. We used th e spec- 
trum of J1148+5251 obtained bv llwamuro et al.l ([2004) 
as our model spectru m. We used the HST Focus Mode0 
(|Cox fc Niemil |2011[ ) to estimate the secondary mir- 
ror despace for each exposure, and included the field- 
dependent coma and astigmatism measurements built 
into TinyTim. The individual HST detectors have differ- 
ent mean focus offsets in the Focus Model, but the offset 
for the WFC3 IR channel has not been characterized. 
We therefore allowed the mean Z4 Zernik e coefficient in 
TinyTim (R® i n the original formulation of lZern ike 1934) 
to float as a free parameter in our optimization, adding 
this single mean value to the Focus Model estimate for 
each exposure. These models for individual exposures 
were then combined, weighted by exposure time, to pro- 
duce a composite PSF for each Multidrizzle-combined 
science image. GalFit also accepts a pixel response con- 
volution kernel for oversampled PSFs. We generated this 
quantity in a similar manner, by drizzling copies of the 



13 http://www.stsci.edu/hst/observatory/focus/TinyTim 

14 http:/ /www. stsci.edu/hst/observatory/focus/FocusModel 
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Fig. 2. Model PSF subtraction. Left panels: Multidrizzle-combined F125W (J, top) and F160W (if, bottom) WFC3 images 
of Jl 148+5251. Middle Panels: TinyTim models of the quasar point source, constructed by optimizing parameters for the PSF star 
observations, then scaled and shifted by GalFit. Right panels: Fit residuals, showing no significant detection of the underlying galaxy 
beyond (X'3 radius. The over- subtracted flux in the central (X'3 (inner circle) occurs because the best-fit model PSFs have more power in 
the central peak than the observations, and is also seen in residuals when modeling the PSF star. This region was excluded from the fit. 
The noise floor in the residual panel is 40% that of the residual panel in Figure ^ From r = (X'3— C/5 (between inner and outer circles) we 
measure a limiting surface brightness of /ij > 23.5, hh > 23.7 mag arcsec -2 (2 a). The noise in the quasar image and uncertainties in the 
PSF model contribute roughly equally within this region. Visible residuals (diffraction spikes and spots) are also seen when subtracting 
this model from the PSF star image. 



empirical WFC3 pixel response convolution kernel (mod- 
eling inter-pixel capacitance and jitter, see Hartig 2008) 
using the same shifts applied to the real images. 

The result of the TinyTim PSF subtraction, with a 
significantly reduced noise floor compared to the direct 
subtraction, is shown in Figure O No host galaxy is de- 
tected, to a limiting surface brightness from r = 0'/3 to 
0'/5 radius of \ij > 23.5, \±h > 23.7 mag arcsec -2 (2 a). 
The inner 0'/3 was excluded from the fit, as the best-fit 
TinyTim models produce PSF cores that are consistently 
narrower than those observed. Visible residuals (diffrac- 
tion spikes and spots) are also seen when subtracting this 
model from the PSF star observations. 

4. HOST GALAXY SIMULATIONS 

Having established no host galaxy detection using the 
TinyTim PSF, we sought to quantify this subtraction 
method's ability to recover the host galaxy flux as a func- 
tion of host galaxy parameters. To do so, we used GalFit 
to simulate a point source along with a Sersic profile host 
galaxy, with total flux adding up to mj = 19.1 mag, 
the measured flux in F125W. This simulated image con- 
tained no noise, so we added both shot noise from the 
object and a Gaussian noise field drizzled in the same 
manner as the real quasar image, to match its correlated 
noise properties. We then performed the same analysis 
that we used on the real quasar image, using GalFit to 
subtract a Tiny Tim-generated point source and measur- 
ing the surface brightness from r = 0'/3 to 0'/5 in the 
residual image. 

We ran a grid of 256 models using this technique, vary- 
ing the total integrated flux of the host galaxy from 



mj = 20—26 mag, the effective radius from r e = O'/l— Of! 9, 
and Sersic indexes n = 1.0 and 4.0. The magnitude range 
represents host galaxies with luminosities from ~ 1/2 to 
1/500 of the total quasar luminosity. Fainter host galax- 
ies than this are undetectable due to shot noise from the 
point source. The range in effective radius corresponds 
to r e ~ 0.6-5.0 kpc at z = 6A2. 

Figure [3] summarizes the results of these simulations, 
plotting the measured surface brightness from r = Of! 3 — 
0% and contours representing the 1,2, and 5 a detection 
limits. Inspecting the residuals of these model subtrac- 
tions, we also found that bright (mj < 22.5 mag), com- 
pact (r e < Of! 3) host galaxies would cause the method 
to significantly over-subtract the PSF. This would show 
negative residuals from the diffraction spikes, which are 
not seen in Figure El 

The model surface brightness from r = 0f!3-0f!b reaches 
the 2 a upper limit of jij > 23.5 mag arcsec -2 for a host 
galaxy of m j > 22—23 mag, depending upon Sersic index 
and effective radius. 

5. DISCUSSION 

We have performed point source subtraction on the 
z = 6.42 quasar Jl 148+5251, with both empirical and 
modeled PSFs. Using direct subtraction, we measure 
an upper limit of mj > 22.8 mag, run > 23.0 mag 
(2 a). With the modeled PSF subtraction we measure 
a limiting surface brightness measured from Of! 3 — Of! 5 of 
jij > 23.5 mag arcsec -2 , \±h > 23.7 mag arcsec -2 (2 a). 
Performing the same subtraction method on simulated 
quasars, we found that this surface brightness limit cor- 
responds to a host galaxy of m j > 22—23 mag, consistent 
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Fig. 3. — PSF subtraction on simulated host galaxies. Surface 
brightness predicted from r = CK'3— (X'5 as a function of simulated 
host galaxy parameters. The simulated observation is approxi- 
mated as a PSF component along with a Sersic profile, with total 
flux adding up to the observed mj from our WFC3 image. The 
integrated magnitude (mh ost ), effective radius (r e ), and Sersic in- 
dex (n = 1.0, left panels and n = 4.0, right panels) of the Sersic 
profile are varied with each model. The 1, 2, and 5cr detection 
significance levels are plotted as black lines. The measured surface 
brightness reaches the 2 a limit (fij > 23.5 mag arcsec - 2 ) for a 
host galaxy of mj > 22 — 23 mag, depending upon Sersic index 
and effective radius. 



with the direct subtraction limit. 

Using the direct subtraction limits, the upper limits 
on the rest-frame monochromatic luminosity (XL\) at 
1700 A and 2200 A are Li 700 < 8.4 x 10 11 L© and L 22 oo < 
5.4xlO n L , assuming a flat spe ctrum within each ban d 
when applying the K-correction (|Oke fc Sa ndage 1968). 
This is comparable to the most lumino us Lyman break 
galaxies at z~2 — 3 (jHoopes et al.ll200"7l ). 

Usin g our upper limits and Equation 1 from iKennicuttl 
(1998), which relates L v to star formation rate, we esti- 
mate a star formation rate of SFR < 210 — 250 M© yr _1 . 
This estimate ignores dust attenuation and assumes 
a continuous star formation rate over 10 8 years or 
longer. A younger population would decrease this upper 
limit, while dust would allow for a higher (absorption- 
corrected) rate. The star formation r ate estimated from 
the AGN-corrected FIR luminosity bv lWang etHI (|2010f ) 
is 2380 M yr- 1 . Since J1148+5251 would be classified 
as a ULIRG locally, this discrepancy is likely due to sig- 
nificant UV absorption by dust. 

We can constrain the infrared excess (IRX) of the host 
galaxy, defined as the infrared to f ar- ultraviolet (FUV ) 
luminosity ratio Lir/Lfuv (e.g., iHowell et al.l 120101 ), 
usually expressed in logarithmic units. Using our upper 
limit for L1700 and an A GN-corrected infr ared luminos- 
ity L IR = 9.2 x 10 12 L© (jWang et al.ll2010D . we calculate 
\og(IRX) > 1.0, consistent with local luminous inf rared 
galaxies (LIRGs) and ULIRGs ([Howell et all 12010( 1. but 
greater than local sta rburst galaxies and high-redshift 
Lyman break galaxies (|Overzier et al.ll2011l ). 

Figure |4] plots broad-band measurements for 
Jl 148+5251, as well as our upper limits for the 
host galaxy flux at 1700 A and 22 00 A and th e AGN - 
corrected FIR measurements of I Wang et"aTl (|2010h . 
Also plotted are the spectral energy distributions 
(SEDs) of four local galaxy systems — three LIRGs 
(Arp 220, IRAS 22491-1 808, and IC 1623), re presenting 
the range in IRX from IHowell et all (|2010f ). and the 
star- forming spiral NGC 4631, representing a galaxy 
with log(IRX) < 1.0, which is thus excluded as a poten- 
tial host. Photometric points for the local galaxies are 
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Fig. 4. — Comparison of broad-band photometry for J1148+5251 
to local galaxies. Blue circles are broad- band (quasar an d 
host) photome t ry of J 1148+5251 t aken f rom Fan et al. (2003), 
l lwamuro ei~aT] (l200l ). I Jiang et all (120 06). Beelen et al. (2006), 
Robson et al. (2004), and Bertoldi et al. (2003). Red squares show 
our upper limits for the host galaxy flux at 170 A and 2200 A, an d 
the AGN-corrected FIR measurements from Wang et al. (2010). 
The light gray spectr u m is the average radio-quiet quasar spec- 
trum of Shan g et al.l (|2011l ), normalized to J1148+5251 from 
0.1 - lfim. The dotted purple SED is NGC 4631, a local spi- 
ral with \og(IRX) < 1.0. Other SEDs are those of the local 
LIRGs Arp 220 (solid red, log(IRX) = 3.423), IRAS 22491- 
1808 (dashed green, log(IRX) = 2.198), and IC 1623 (dot-dashed 
blue, log(IRX) = 1.3 79), representing h igh, average, and low IRX 
LIRGs, respectively (Ho well et alJ feOlO). The local galaxy SEDs 
have been normalized to match the AGN-corrected emission of 
J1148+5251 between rest-frame 40 anc! 200 microns. Our con- 
straint of \og(IRX) > 1.0 (and the 2200 A flux limit) matches most 
local LIRGs, but is greater than almost all local star- forming galax- 
ies and high-r edshift Lyman break galaxies (Howell et al. 120101 . 
IQverzier et aI1[201ll ). 

taken from NEEF^l and the SEDs have been normalized 
to match the AGN-corrected emission of Jl 148+5251 
between 40 and 200 microns. 

We estimate Apuv > 2.1 mag of UV absorption using 
the relation between the IRX flux ratio and Apuv (e.g., 
IQverzier et al. I [20 111 Equation 1). Using the empirical 
relation IRX M m.i nn ^ r (A FUV = 4.54 + 2.07/3 + 0.4) from 
IQverzier et al.l (|2011l ) , we obtain a limit of > — 1.2 + 
0.2. This matches local (U) LIRGs ([Howell et al.ll2010D . 
but is redder than almost all local star-forming; galaxies 
and z ~ 6 Lyman b reak galaxies (jOverzier "etatl[20TTL 
iBouwens etal]l2Q12f l. 

The Tiny Tim-based subtraction may be improved in 
the future with more accurate WFC3 IR PSFs. Since 
uncertainties introduced by the PSF model scale with 
PSF brightness, our further WFC3 observations will tar- 
get quasars where the contrast ratio between point source 
and host galaxy is expected to be smaller, such as opti- 
cally faint z~6 quasars with large FIR luminosities. The 
James Webb Space Telescope will enable us to use the 
PSF subtraction method at rest-frame ultraviolet and 
optical wavelengths with better-sampled empirical PSFs 
in a more stable thermal environment. 

Special thanks are due to K. Long, J. MacKenty, and 
A. Roman at STScI for their advice and expertise when 

15 The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Technol- 
ogy, under contract with NASA. 
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planning the observations, and the anonymous referee 
for very useful suggestions. We are grateful to J. Krist 
and C. Peng for the public availability of Tiny Tim and 
GalFit, respectively. This research has made use of 
the SIMBAD database, operated at CDS, Strasbourg, 
France. Support for this HST program was provided 



by NASA through grant GO-12332.*.A from the Space 
Telescope Science Institute, which is operated by the As- 
sociation of Universities for Research in Astronomy, Inc., 
under NASA contract NAS 5-26555. We owe a deep debt 
of gratitude to the entire WFC3 instrument team and the 
astronauts of STS-125 for this marvelous instrument on 
HST. 
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